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A solution has been obtained for the complete tunnel- interference 
flow for a lifting vortex in a two-dimensional slotted tunnel. Curves 
are presented for the longitudinal distribution of tunnel- induced down- 
wash angle for various values of the boundary openness parameter and for 
various heights of the vortex above the tunnel center line. Some quanti- 
tative discussion is given of the use of these results in calculating 
the tunnel interference for three-dimensional wings in rectangular tun- 
nels with closed side walls and slotted top and bottom. 


INTRODUCTION 


The problem of determining lift interference corrections for slotted- 
throat tunnels by using an approximate "homogeneous" boundary condition 
has been treated by several authors (refs. 1 to 5). For both two- and 
three-dimensional tunnels, results have been given for the vertical inter- 
ference velocity at a lifting vortex located on the tunnel center line. 
(See, for example, ref. 1.) 

This report extends the two-dimensional results of reference 1 by 
giving the longitudinal distribution of this vertical interference veloc- 
ity along the horizontal line through the vortex, for various vertical 
positions of the vortex in the tunnel. The manifest application of these 
results is to the determination of the tunnel- induced angle of attack 
and camber in two-dimensional airfoil testing. An important further appli- 
cation is to the determination of the tunnel-induced downwash angle in 
three-dimensional-model testing. The feasibility of the latter appli- 
cation, was pointed out in reference 2 and will be further discussed in 
the present report. 

As in the previous tunnel-interference studies, the analysis is made 
for incompressible flow. For the present problem, which considers only 
the tunnel- induced vertical velocity for given lift, the compressibility 
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ef fect can be taken into account merely by reading the induced angle at 
x/l 1 - M 2 instead of at x (Gothert rule) . 

The method of solution of the problem is considered to be of partic- 
ular 1 interest. It uses a generalized Fourier series in a form that 
automatically satisfies the boundary conditions. The method has appar- 
ently not been used in previous wind-tunnel . interference studies, although 
similar methods of solution have been used, for example, for problems 
involving thermal diffusion (ref. 4 ) and stress diffusion (ref. 5). 


SYMBOLS 


^ constant coefficients in the series expansion for cp 



c chord of two-dimensional airfoil 

lift coefficient of two-dimensional airfoil 

C area of tunnel cross section 

Cl wing lift coefficient 

d slot spacing (distance between center lines of slots) 

k height of vortex above tunnel center line, expressed as 

fraction of tunnel semiheight 

h tunnel semiheight 

l apparent mass associated with flow normal to slotted 

boundary, per unit area of boundary, ^ log esc £2. 

’ it e 2 

M Mach number 

r n root of equation (2a) 

R n root of equation (2b) 


S 


wing area 
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V tunnel free-stream velocity 

x, z rectangular coordinates, expressed as fraction of tunnel 

■semiheight (fig. l) 

T vortex strength, positive clockwise (lifting) 

5 tunnel correction factor 

^2-dim two-dimensional- tunnel correction factor 

8 C that part of the total tunnel correction factor contributed 

by the center row of images 

AB error in tunnel correction factor incurred by using two- 

dimensional vortices in place of the exact images in the 
rows above and below the center row; defined only for 
closed tunnels or for tunnels with completely open top 
and bottom boundaries 

e tunnel- induced downwash angle 

9 total perturbation potential (due to lifting vortex and 

tunnel interference) 

9p potential of lifting vortex 


9 2 potential of tunnel interference 

a open ratio of slotted boundary (ratio of slot width to 

slot spacing) 


BASIS OF METHOD OF SOLUTION 

The boundary condition at a slotted wall, considered as an equiva- 
lent homogeneous boundary, is derived in references 1, 2, and 5 as 
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perturbation potential 

derivative of cp in the direction of the outer normal 
at the boundary 

l = — logg esc — 

it e 2 

d slot spacing 

a open ratio of slotted boundary (ratio of slot width to 

slot spacing) 

For the present problem of the two-dimensional tunnel it will be 
convenient to make all lengths nondimensional by dividing by the tunnel 
semiheight h; in other words , the tunnel semiheight will be taken as 
unity. A sketch of this tunnel showing the coordinate system and the 
lifting (bound) vortex is shown in figure 1. The boundary condition 
may be written 

9 


9 

where g = 

h 

The boundary condition can be automatically satisfied by making 
use of a trigonometric Fourier series of a more general form than that 
commonly used. Consider the functions sin r n z and cos I^z 

(n = 1, 2, J>, . . .) where the r n are the positive roots of 

tan r + gr = 0 (2a) 

and the R n are the positive roots of 


+ g = 0 (z - 1) (la) 

oz 

- 8 I 2 » 0 (z = -1) (lb) 

oz 


4 

where 

9 

8cp 

8n 
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The functions sin r n z and cos R n z form a complete orthogonal set 
over the interval -1 S z = 1 (ref. 4). The series 


s 


00 

Z An sin r n z + Bn cos R n z 
n=l 


satisfies the conditions (la) and (lb) at z = 1 and z = — 1, respec- 
tively, because 


s + 


g ^ 
8z 


Z=1 


00 

= ^ V^sin r n + gr n cos r n ) + B^cos Rn - gR n sin Rn) 


n=l 



1 


00 

^ A n (-sin r n - gr n cos r n ) + Bn (cos R n - gR n Bin R n) 
n=l 


where the quantities in parentheses are all zero because of conditions (2a) 
and (2b) . Hence, any series of the form 


s 


00 

^ A^^sin r n z + B^^cos R n z 
n=l 


representing a function in the xz -plane automatically satisfies the 
boundary conditions (la) and (lb) . 

It is of interest to note that the arguments of the functions of a 
Fourier -Bessel series can be similarly chosen so that each term satisfies 

the boundary condition cp + g = 0 on the boundary of a circle (ref. 6) . 

5n 

APPLICATION TO THE PROBLEM OF THE LIFTING VORTEX 


Bonriflary conditions .- In reference 1 it was shown that the total 
perturbation flow (vortex flow plus boundary interference) in the tunnel 
could be considered to have two parts: (l) an antisymmetric part 
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characterized "by a uniform upwash velocity far upstream, an equal uniform 
downwash velocity far downstream, and zero downwash along the vertical 
axis (the z-axis) through the vortex itself, and ( 2 ) a flow with uniform 
downwash velocity throughout the entire tunnel, of such magnitude that 
it cancels the upstream upwash velocity of part (l) and thereby provides 
zero net effect far upstream and uniform downwash velocity ■ along the 
z-axis. The value of this uniform downwash velocity (for the tunnel of 

unit semiheight) was shown to be — . This value was derived in 

4(g + 1 ) 

reference 1 for the case of the vortex at the origin; however, the line 
of reasoning is not altered by changing the vertical location of the 
vortex in the tunnel, so that the same expression remains applicable in 
the general off center case. Accordingly the potential along the z-axis is 

given by £2 f £ above the vortex location and by - Ez £ 

^(g +1) 4 4(g +1) 4 

below the vortex location. (The potential is here considered to increase 
in the direction of flow, contrary to the usage in reference 1; also, a 
clockwise (lifting) vortex is here considered as positive.) 

The problem will be solved for the upstream half of the tunnel: 
x ^ 0, -1 S z = 1; the solution for the downstream half follows from 

the known symmetry characteristics. For the upstream half, then, the 
problem, as given by the preceding discussion, is to find an expression 
for 9 that satisfies Laplace's equation, satisfies equations (la) and 
(lb), approaches zero as x approaches and along the z-axis equals 

£2 + £ above the vortex and £2 £ below the vortex. 

4(g +l) 14 4 (g +1) 4 


Solution for total perturbation flow .- The total perturbation 
potential is expressed in the form 


r n x 


RnX 


cp = ^ A n e “ sin r n z + B^e “ cos R n z 
n=l 


(3) 


where the r n and R n are the positive roots of equations (2a) and (2b), 
respectively. The coefficients A n and 3^ are determined by matching 
the expression with the known potential along the z-axis . That is, let 

r 

0 


A n sin r n z + Bn cos R n z = - 


Fz 


n=i 


4(g + 1) 




r 


(k < z < 1) 
(z = k) 
(-1 < z < k) 
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Since the terms of the series are orthogonal in the range (-1, l ) , each 
A n is found by multiplying both sides of the equation by sin r n z dz 

and integrating from -1 to 1: 



r 

Ms + 1) 
r 

Mg + 1) 



p 

sin r n z dz - — 
k 

sin r n z dz + E 



sin r n z dz + 



sin r n z dz 



sin r n z dz 


Applying the same procedure to the cosine term gives 




cos^RqZ dz = - 


r 

£ 


p k p 1 

J cos RjjZ dz + E J cos B^z dz 


r \ 

f 

2 J 0 


cos R n z dz 


Performing the indicated integrations and simplifying the results by 
means of equations (2a) and (2b) yields 


An 



cos r n k 
+ g cos 2 r n 


p sin R n k 
SR ^- 1 + g sin% n 


Substituting these expressions for the coefficients in equation ( 3 ) 
completely determines cp in the region x < 0. 

Tunnel-interference velocity . - The potential cp just derived is 

the sum of the potential cp-i = - — tan”-*- 2 ~ of the lifting vortex 

2rt x 

and the potential cp 2 of the tunnel interference: cp = cp-j_ + cp 2 . The 

dcpo 

desired tunnel-interference upwash velocity is found by subtracting 

oz 
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^2 8cp ^q>i 

3z 8z 8z 


f r n x R n x \ 

e cos r n k cos r n z ^ e sin R n k sin R n z \ 

V 1 + g cos 2 r n 1 + g sin 2 R n / 


(x < 0) 


The entire model may generally be considered with adequate accuracy to 
lie in one horizontal plane. Then, since the interference is of inter- 
est only in this plane, z may be set equal to k in the preceding 
equation to give 


r x 

2jt X 2 ■+ (z - k) 2 



&p 2 r y / e rnX cos 2 r n k + e RnX sin 2 R r ik \ + T 

Z 2 ^ l 1 + g cos 2 r n 1 + g sin ^R n J 2nx 


(x < 0) (4) 


Equation (14-) gives the tunnel-induced upwash velocity in the 
upstream part of the tunnel. At x = 0, the upwash velocity is 

p 

- — as already mentioned. The tunnel-induced upwash velocity 

Mg + 1) 

in the downstream part of the tunnel follows from the symmetry charac- 
teristics previously described; that is. 


M 

1 

u 

= ^2 


r 

4(g + 1 ) 8z J 

Sz 

X, z 

-X, z 

L Mg + 1)J 


or 



r 

2(g + 1) 



z 


where the symbol 
indicated by the 




means the value of 
subscript. 


the derivative at the point 
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Calculated results .- The tunnel-induced vertical velocities were 
computed for points along the horizontal line through the hound vortex, 
for ranges of values of g and k that are considered to more than 
cover the likely ranges of slot designs and airfoil locations. The 

results are plotted in figure 2 in nondimensional form as where 

e is downwash angle and V is tunnel speed. For direct application 
to two-.dimensIonal airfoil testing, the equivalent expression in terms 

of chord c and section lift coefficient c 7 , is also shown. The 

cc y 

expression corresponds to the usual three -dimensional- tunnel correction 
factor 8 defined by 


e = 8 | C L 


where S is wing area, C is the area of the tunnel cross section, and 
is the wing lift coefficient. 

The two limiting cases - the closed tunnel (g = °o) and the open 
tunnel (g = 0) - are not included in figure 2 because closed expressions, 
readily derived by the method of images, are available for these cases. 
For the closed tunnel. 


8 


2 -dim = 


8 


cosh cos 2 Sjl 
4 2 


sinh ( B inh 2 £2L + cos 2 — 

k V 4 2 / it- J 


For the open tunnel. 


ft _ 1 . 1 

°2-dim “ 4 + 8 


cosh 




itx 

k 


ESTIMATING CORRECTIONS FOR THREE-DIMENSIONAL MODELS 


Suggested correction as the sum of two parts .- In reference 2 it 
was suggested that two-dimensional results such as those derived in the 
preceding section could be applied in determining corrections for three- 
dimensional models in rectangular wind tunnels with closed side walls 
and slotted upper and lower boundaries, provided the height-width ratio 
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of the tunnel Is not too small. Specif ically, the tunnel-induced angle 
at any point would he the sum of two parts: (l) that due to the infinite 

row of images of the model in the vertical walls, and (2) the two- 
dimensional corrections of figure 2, determined as for a two-dimensional 
wing of the same lift at the same longitudinal location in the tunnel. 

The purpose of the present section is to review the theory of this method 
and to provide some quantitative discussion of the method and of the tun- 
nel height-width ratios for which it can he considered satisfactory. For 
the smaller height-width ratios, for which the method is not quite valid, 
a further suggestion is given hy which the method might still he used 
with an additional modification. 

A complete analysis should include comparisons between induced angles 
determined hy this suggested method and those determined hy some exact 
method similar, say, to those outlined in reference 1. Such exact three- 
dimensional calculations would he very laborious, however, for the caBe 
of the tunnel with slotted upper and lower boundaries. Accordingly, the 
present approach will he to make exact calculations for the cases of com- 
pletely open and completely closed top and bottom boundaries, which are 
the two limiting cases of slotted boundaries, and then to apply arguments 
of plausibility in extending the results to the case of slotted boundaries. 

Formulation of problem for the two limiting cases .- Figure 3(a) 
shows the . cross section of a closed rectangular tunnel (heavy outline) 
with its nearest images. The plus sign in the center of the tunnel 
represents the semi-infinite doublet line (or horseshoe vortex of zero 
span) trailing downstream from a point concentration of lift in the cen- 
ter of the tunnel. The image doublet lines are indicated by plus or 
minus signs according as they are the same as or the reverse of the cen- 
ter one. The tunnel-induced downwash or upwash velocity at any point in 
the tunnel (say at the lifting element itself or at some downstream point 
representing the tail location) is the net induced velocity due to all 
the images. The concepts under lying this figure and the application of 
this theory have been widely used in wind- tunnel-correction analysis (for 
example, refs. 7 and 8) and will not be further discussed here. 

As pointed out in references 7 and 8, the field of a doublet line 
at a point sufficiently distant from it is practically the same as the 
field of a narrow horseshoe vortex having the same moment (same product 
of vortex strength and span - in other words, the same lift) . Suppose, 
then, it were valid to replace all the image doublet lines in the rows 
above and below the center row by horseshoe vortices of span equal to 
the tunnel width. The trailing vortices would all cancel in pairs and 
only the bound vortices, now joined to form continuous infinitely long 
vortices, would remain, as indicated in figure 3(t>) • This set of joined 
bound vortices, however, will be recognized as identical with the set of 
image bound vortices that are used to calculate corrections for a closed 
two-dimensional tunnel (ref. 9) • If the image system of figure 3("b) were 
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indeed an accurate approximation for the image system of figure 3(a ) , 
it would follow that the tunnel-induced angle for this case is the sum 
of (l) the effect of the images in the center row and ( 2 ) the two- 
dimensional correction, just as suggested in the opening paragraph of 
the preceding section. 

The corresponding image diagrams for the tunnel with open top and 
bottom are shown in figures 4(a) and Mb). The argument is essentially 
the same as for the closed tunnel and will not be repeated. 

The problem now is to determine for what height-width ratios the 
image systems of figures 3(b) and Mb) satisfactorily approximate the 
exact image systems of figures 3(a) and 4(a), respectively, so that the 
suggested calculation procedure would be accurate. 

Calculations of the error .- In order to investigate this problem, 
calculations were made of the error in the induced-angle factor S along 
the tunnel center line that would result from using the image systems of 
figures 3(b) and 4(b) instead of those of figures 3(a) and 4(a), respec- 
tively. The calculations were made for tunnel height-width ratios of 0.6, 
0.75, 1.0, and 1.2. These calculations for the zero-span horseshoe vor- 
tices (or point concentrations of lift in the center of the tunnel), how- 
ever, represent the least favorable case. In order to investigate the 
problem for conditions more nearly representative of typical wing tests, 
calculations were also made for horseshoe vortices of span equal to 0.4 
the tunnel width. As a somewhat extreme case, some calculations were 
also made for horseshoe vortices of span equal to 0.8 the tunnel width. 
For the l i m iting case, for which the horseshoe vortex has a span equal 
to the tunnel width, the exact image system would be identical with the 
two-dimensional image system, and thus the error would be zero for any 
height-width ratio. 

The calculated errors are plotted in figure 5 as the error Zf> in 
the usual tunnel correction factor 6 in the equation 

e = 6 |°L 


Specifically, is the induced downwash factor given by the exact 

image system minus the induced downwash factor given by the approximate 
image system. The approximate image system gives too small a downwash 
angle (too small a factor) for the open top and bottom and too large a 
downwash angle (too large a factor) for the closed tunnel. In other 
words, is positive for the tunnel with open top and bottom, and is 

negative for the closed tunnel. The value of £5 at x = 0 represents 
the error at the lifting element itself; values at positive values of x 
represent the error in downwash behind the wing. The upstream (negative x) 
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parts of the curves are not shown, since the curves all have the same 
type of symmetry as the curves of figure 2. 

As an indication of the order of magnitude of the error that might 
arise, consider the largest value of dB in figure 5, 0.092, for the 
tunnel with a height-width ratio of 0.6. If the wing area S is one- 
tenth the tunnel cross-sectional area C and the lift coefficient 

is 1.0; the corresponding error in tunnel-induced angle is 
0.092 x 0.1 x 1.0 x 57.3° = 0.53°. The error decreases rapidly with 
increasing height-width ratio and with increasing span of the horseshoe 
vortex, and it is well within usual experimental accuracy for height- 
width ratios of about 1.0. 

Application to the slotted tunnel .- For those tunnel height-width 
ratios for which the results show no practically significant differences 
between the field of the true image system and the field of the system 
of infinite line vortices that approximates it, it would seem reasonable 
to suppose that the effect of slotted top and bottom boundaries would 
be similarly two-dimensional. In those cases in which the difference 
does become practically significant (for example, for tunnel height- 
width ratios less than 1.0 with the larger lift coefficients and wing 
areas) the matter is less obvious. A suggested approach is as follows: 

Let the curves of figure 5 be considered as correction curves rather 
than as error curves. That is, suppose that the tunnel-induced angle 
(for the completely open or completely closed top and bottom) is con- 
sidered to be the sum of the following three components : 

(1) That due to the images in the center row 

(2) The two-dimensional induced angles corresponding to the infinite 
line vortices of figure 3("b) or *0>) 

(3) The correction indicated by figure 5, which, essentially, cor- 
rects the two-dimensional values of component (2) to the correct three- 
dimensional value corresponding to figure 3(a) or 4 (a). 

Then a parallel procedure for the tunnel with slotted top and bottom 
would be to suppose the tunnel-induced angle to be the sum of the fol- 
lowing three components: 

(1) That due to the images in the center row 

(2) The two-dimensional induced angles given by figure 2 for the 
appropriate value of g 

(3) A correction determined by interpolating between the closed 
and open cases of figure 5. The interpolation could be accor ding to 
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the "degree of openness" of the slotted "boundary as indicated, say, by 

the factor — - — that determines the two-dimensional correction at 
B + 1 

x = 0. This factor is zero for the closed tunnel (g = ») and 1 for the 
open tunnel (g = 0) . According to this suggested interpolation proce- 
dure, the contribution of this third component to 8 is 

^closed + g +. 1 (^®open “ ^closed) 

where closed and - ^open are the values of Zfi given by figure 5 
for the closed and open cases, respectively. Furthermore, since, as 
figure 5 shows, ^B open and closed are a] - mos ' fc equal, except for 

sign, this expression may be approximated as 


open g 4 - 1 open 


or 


1 + g open 

If the slot design is such that g = 1, this third component will be 
approximately zero. 

In order to help in the application of thiB procedure, the first 
component 8 C , due to the images in the center row, has been plotted in 
figure 6 . The net induced downwash factor, finally, is 

8 = 5 C + S 2 _ dim + i-=-| £B open 

where S c is given by figure 6 , 8 2 _ dim is given by figure 2 , and £ 5 open 
is given by figure 5 (positive values). 

Asymmetric locations of the model in the tunnel .- A sting-mounted 
model might, at high angles of attack, be considerably above the t unne l 
center line. Figure 7(a) shows the image system for such an asy mme tric 
case, for a point concentration of lift in a closed t unn el - The corre- 
sponding approximation by means of infinite vortices is shown in fig- 
ure 7(b) . No calculations were made for this case, but it is apparent 
from the image system that the error is very nearly the average of t ha t 
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for a t uim el whose height is the distance m-^ to the nearer row and that 

for a tunnel whose height is the distance to the farther row. Except 

for this one modification, the discussion for the slotted-boundary tunnel 

follows as before. 

Yawing a sting -mounted model also tends to increase the error. Fig- 
ure 8(a) shows the Ima ge system for an extreme case in which the point 
concentration of lift has been brought close to one of the side walls. 

The image system in this case is very nearly the same as for a centered 
model in a tunnel of half the height-width ratio (fig. 8(b)), for which 
the error is greatly increased. Fortunately, lateral displacement of 
the lifting point to the extent shown in figure 8(a) is far in excess 
of what might be considered a realistic displacement due to yawing a 
model in a wind tunnel. 


CONCLUDING- REMARKS 


A practical mathematical method is described for determining the 
flow in a two-dimensional slotted wind tunnel containing a lifting wing. 
Calculated results are given for the longitudinal distribution of tunnel- 
induced angles for various slot geometries and for various vertical loca- 
tions of the wing in the tunnel. Some quantitative discussion is given 
of the application of these results to the testing of three-dimensional 
lifting models in a tunnel with closed side walls and slotted top and 
bottom. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 6, 1958- 
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Figure 2.- Longitudinal distribution of tunnel- induced downwash angle 
for a lifting vortex in a two-dimensional slotted tunnel. 
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Figure 2.- Continued 
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(d) — - — = 0.6. 
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Figure 2.- Continued. 
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Figure 2.~ Continued 
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Figure 2.- Continued 
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(h) 
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Figure 2.- Continued. 
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Longitudinal distance from vortex, tunnel semiheights 

(i) — - — = 0.1. 

S + 1 


Figure 2.- Concluded 
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(a) Exact Image system. 



(b) Approximate image system. Doublet lines above and below the center 
have been replaced by horseshoe vortices of span equal to the tun- 
nel width. 

Figure 3.- Image systems for a doublet line in a closed 

rectangular tunnel. 
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(b) Approx: 
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Figure 4.- 
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Longitudinal distance from lifting element, tunnel widths 


Figure 5.- Plots of £8 (correction factor given by the exact image 
system of fig. 3(a) or 4(a) minus correction factor given by the 
approximate image system of fig. 3(b) or 4(b)). Positive values are 
for the tunnel with open top and bottoml and negative values are for 
the closed tunnel. 
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(a) Doublet near the side wall. 
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(b) Doublet at center of a tunnel having twice the width of the tunnel 

diagrammed above . 

Figure 8.- Image system, for a doublet near the side wall of a closed 
rect ang ular tunnel, and comparable image system for a doublet at the 
center of a tunnel twice as wide as the first. 
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